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Abstract 

We develop the perturbative QCD formalism including Sudakov effects 
for semi-leptonic B meson decays. We evaluate the differential decay rate of 
B — > ttIv, and find that the perturbative calculation is reliable for the energy 
fraction of the pion above 0.3. Combining predictions from the soft pion 
theorems, we extract the value of the matrix element \V u b\ which is roughly 
2.7 x 10~ 3 . 



Exclusive semi-leptonic meson decays have been intensively studied by 
several approaches, which provide the information of the mixing angles in the 
Cabibbo-Kobayashi-Maskawa matrix of the standard model. Chiral symme- 
try and heavy quark symmetry have been applied to the decays K — > tiIu [|TJ 
and B — > Dlv p|, respectively, from which model- independent extraction of 
the matrix elements \y us \ and |T4&| is obtained. For the heavy-to-light transi- 
tion B — > nlu, which gives reliable estimation of \ V u b\ ||, neither of the above 
theories is appropriate. Recently, this decay has been investigated using the 
heavy quark effective theory (HQET) || , in which only the normalization of 
the relevant form factors is determined in terms of the soft pion relations. 
An explicit evaluation of the decay B — > nlu based on the perturbative QCD 
(PQCD) formalism has been proposed ||. However, it leads to results which 
are too small compared to current experimental data. 

In this letter we shall develop a modified PQCD approach to heavy me- 
son decays, which includes Sudakov effects || 0]. These effects, arising from 
the all-order summation of large radiative corrections, suppress contributions 
from the nonperturbative region, and have been found to extend the appli- 
cability of PQCD down to the energy scale of few GeV in the study of elastic 
hadron form factors 0. We shall show that the PQCD approach is proper 
for the decay B — > -kIu, at least when the pion is energetic. 

The amplitude of the considered process is written as 

A(P U P 2 ) = ^V^{1 - l5 )l{i:{P 2 )\u^b\B{P 1 )) , (1) 

where Gp is the Fermi coupling constant, and Pi and P 2 are the B meson and 
pion momenta, respectively. We start with the lowest-order factorization for 
the matrix element M M = (it (P 2 )\u / ~f fJ 'b\B(Pi)) , in which the b quark carries 
the momentum P\ — k\, and its light partner carries h\. These momenta 
satisfy the on-shell conditions (Pi — ki) 2 = m 2 , Pj 2 = m 2 B and k\ = 0, 
and tub being the b quark and B meson masses, respectively. We choose 
the Breit frame such that P x + = Pf = m^/-\/2, P 2 = r]niB/V2 and all 
other components of Pj's vanish, where r] is related to the energy fraction of 
the pion by P° = rjmB/2, < r] < 1. ki has a minus component, defining 
the momentum fraction X\ = fcj~/Pf in the B meson wave function, and 
small amount of transverse components kix- The light valence quark of 
the B meson, after absorbing the hard gluon, goes into the pion with the 
momentum fraction x 2 and transverse momenta h 2 T- 
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We then consider how to group radiative corrections into the basic factor- 
ization by locating their leading momentum regions, from which important 
contributions to loop integrals arise. The important corrections are charac- 
terized by large single logarithms, which are either collinear or soft. These 
two regions may overlap and give double logarithms. It is known that single 
logarithms can be summed to all orders using renormalization group (RG) 
methods, while double logarithms must be organized by the resummation 
technique ||, which has been developed in axial gauge n ■ A = 0, n being the 
gauge vector and A the gauge field. 

A careful analysis shows that reducible corrections on the pion side pro- 
duce double logarithms with soft ones cancelled in the asymptotic region 
b — > j7|, b being the conjugate variable to ky. Hence, they are dominated 
by collinear enhancements, and can be absorbed into the pion wave funtion. 
Reducible corrections on the B meson side also give double logarithms, but 
the soft ones do not cancel and the collinear ones are suppressed by the B me- 
son wave function. Therefore, these corrections can be absorbed into the B 
meson wave function, which is also dominated by soft dynamics. Irreducible 
corrections, with an extra gluon connecting the pion and the B meson, give 
only soft divergences, which cancel asymptotically. They are then absorbed 
into a hard scattering amplitude. Hence, the factorization picture holds after 
radiative corrections are included. 

With the above reasoning, the factorization formula for M M in the trans- 
verse configuration space can be written as, 

x H> x (x 1 ,X2, bi, b 2 , m, //) V B {xi, bi, Pi, //) , (2) 

with V n and Vb the pion and B meson wave functions, respectively, and 
H** the Fourier transform of the hard scattering amplitude to b space, /i 
is the factorization and renormalization scale. Both V n and Vb contain 
double logarithms, which will be summed up below. The approximation 
mj, m tub = m = 5.28 GeV has been made to simplify the analysis. 

We outline the resummation procedure employed for V w H [jj . If the dou- 
ble logarithms are grouped into an exponential V ~ exp[— In m ln(lnm/ In b)} , 
the problem will become simpler by considering the derivative dP/dlnm = 
CV, where the coefficient C contains only single logarithms, and can be 
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treated by RG methods. Because of the scale invariance of n in the gluon 
propagator, 



N» u (q) 



-i 
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_ n"q" ■ <,",," 2 q*q v 



n ■ q 



(n ■ q)' 



(3) 



V n depends only on the ratio v\ = (P 2 ■ n) 2 /n 2 . It is then possible to relate 
the derivative dP 7r /dlnP 2 + to dV^/dn, which can be easily computed using 
the relations dN^ v /dn a = —(N^ a q u + N ua q^)/q ■ n. The momentum q ap- 
pearing at both ends of the differentiated gluon line hints the application of 
the Ward identity After adding together all diagrams with different differ- 
entiated gluon lines, we obtain an equation graphically described by fig. la, 
in which the square vertex represents gT a v?P^ '/ '{P 2 -nq-n), T a being related 
to the Gell-Mann matrices A a by T a = \ a /2. 

Due to the factor 1/q-n in the new vertex and the nonvanishing of n 2 , the 
leading regions of q are soft and ultraviolet, in which fig. la can be factorized 
according to fig. lb to lowest order of a s . The part on the left-hand side of the 
dashed line is exactly V n , and that on the right-hand side is assigned to the 
coefficient C. We introduce a function K to organize the soft enhancements 
in the first two diagrams of fig. lb, and Q for the ultraviolet divergences of the 
other two diagrams. The soft subtraction in Q is to avoid double counting. 
We then derive the differential equation, 



dlnP 9 



-V 



£?((! - x 2 )v 2 /ij) 



V 



(4) 



where the functions K, and Q have been calculated using RG methods |J. 
Solving eq. (HI), we obtain the solution 



V n (x 2 ,b 2 ,P 2 ,n) = exp 



8 (t,h,v™) 

£,=X2, 1—X2 



Vir{x2, b 2 , fj) . (5) 



The explicit expression of the exponent s(£, b, Q) has been obtained in |]], 
and will not be shown here due to its complexity. 

The function still contains single logarithms from ultraviolet diver- 
gences, which need to be summed using RG methods M. The large-6 behav- 
ior of V n is then written as 



V w = exp 



J2 s(£,,b 2 ,r}m) -2 [ — j q (g(fi)) 
, J 1/62 /i 



(,=X2, 1-X2 



^(x 2 , b 2 



(6) 
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7g = —a s /iT being the quark anomalous dimension in axial gauge. The bi 
dependence in <p n , which corresponds to the intrinsic transverse momentum 
dependence of the pion wave function ||, will be neglected. 

As to Vb, the resummation of the double logarithms is subtler. The self- 
energy correction to the massive b quark, giving only soft single logarithms, 
should be excluded. On the other hand, Vb involves the invariants such as 
Pj 2 , which can not be related to n, so that the technique of replacing d/dm 
by d/dn fails. However, the above difhculities can be removed by applying 
the eikonal approximation to the heavy quark line. In the collinear region 
with the loop momentum q parallel to k\ and in the soft region, the b quark 
line can be replaced by an eikonal line: 

(/>!- M+ A + m)j a (Pi -£4)° 



+ R 



(7) 



(Pi - h + qY - m 2 {P x -k x )-q 

where the remaining part R is less important. The involved physics is that 
a soft gluon or a gluon moving parallelly to k\ can not explore the details of 
the b quark, and its dynamics can be factorized. 

The first difficulty is then resolved, because self-energy diagrams of an 
eikonal line are excluded by definition fllPf . With the scale invariance of 
Pi — k\ as shown in eq. (]?]), which is equivalent to the flavor symmetry in 
HQET, Pi — k\ does not lead to a large scale, and the remaining large scale 
is only fcf . Furthermore, an explicit lowest-order investigation shows that 
Vb depends only on the single ratio v\ = (ki ■ n) 2 /n 2 , and thus d/d/cf can 
be replaced by d/dn now. 

Following the similar procedures to those for the pion, we obtain the 
differential equation, 



= ^C(M + -Q{Pi/ii)\ V B ■ (8) 

It can be shown that the functions K, and Q for the B meson are exactly the 
same as those in eq. (HI). It is then straightforward to derive the solution 



V 



b = exp 



-s(x 1 , 61, m) - 2 / — n/q(g(v)) 



(j)B(xi,h 



At last, the RG analysis of P M gives 



P M (xj,6i,m,/i) = exp 



* d/i 



-4 / -^ lq ( g m 



H^fa, bi, rn, t) , 



(9) 



(10) 
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where t is taken as the largest mass scale associated with the hard gluon, 
t = max (y'xixi^m, 1/62)- Having factorized all the large logarithms 
into the exponents, we can then compute if M to 0(a s ). 

Substituting eqs. (H), (|D|) and JIHD into (|2|), we obtain the factorization 
formula for = fiP\ + f 2 P 2 , where the form factors f\ and f 2 are given 
by 

fi = -l§TTC F m 2 / dxidx 2 / b 1 db 1 b 2 db 2 (p B (x 1 ,b 1 )(j) n (x 2 ) 
Jo Jo 

xx 1 7]h(xi,X2,bi,b 2 ,m)exp[-S(x i ,bi,m)] , (11) 

and 

rl /*oo 

f 2 = lQnCFm 2 / dxidx 2 / bidb 1 b 2 db 2 (j)B(xi,b 1 )(f) 7r (x 2 ) 
Jo Jo 

x [xihfa, x 2 , 61, b 2 , m) + (3 - x 2 r])h(x 2 , x x , b 2 , bx,m)\ 

x exp[-5'(x i ,6 i ,m)] , (12) 

respectively, with 

h(x 1 ,x 2 ,b 1 ,b 2 ,m) = a s (t)K (^/xix 2 r]mb 2 ) 

x[Q{bi - b 2 )K (^/x^mb 1 )I (^/x~[rjmb 2 ) 

+9(b 2 - 6i)Ko(v / ^ m& 2)/o(v / 3 7 i^ m& i)] • ( 13 ) 

Cp is the color factor defined by tr(T a T a ) = N c Cp, N c being the number of 
colors, and K and I are the modified Bessel functions of order zero. The 
complete Sudakov exponent S is given by 



S(xi,bi,m) = s(xi,bi,m) + s(x 2 ,b 2 ,r]m) + s(l - x 2 ,b 2 ,r]m) 
1 



ln lnft/A) | ln ln(t/A) 



ln(l/6iA) ln(l/6 2 A) 



(14) 



with (3 = (33 — 2n/)/12, rif = 4 being the number of quark flavors, and 
A = Aqcd = 100 MeV here. The Sudakov factor exp(— S) decreases quickly 
in the large bi region and vanishes as bi > 1/A. We have kept the intrinsic 
transverse momentum dependence in Since the B meson wave function 
is dominated by soft dynamics, this dependence is more important than that 
in the pion. 
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^ is chosen as the Chernyak-Zhitnitsky model [ OJ , 



<f> v {x) = 5V3f n x(l - x)(l - 2 



x 



(15) 



with f n = 93 MeV the pion decay constant. For the B meson wave function 
we consider |T2| 



C + 



m 



+ 



1 — x x(l — x) 



(16) 



with N = 1.232 GeV 3 and C = -0.99993m 2 , which are obtained from the 
normalizations J dx J d 2 k T $ B = /b/2a/3 and / dx J d 2 k T $| = 1/2, f B = 
160 MeV being the B meson decay constant ||13|| . The Fourier transform of 
<&b gives 



<f> B (x,b) 



irNbx 2 (l 



x) 



m 2 x + Cx{l — x) 



K\ Jm 2 x + Cx{l — x)b 



(17) 



K\ being the modified Bessel function of order one. 

Results of f\ + /2 with b\ and 62 integrated up to the same cutoff b c 
are shown in fig. 2a. We observe that at 77 = 0.3 approximately 50% of 
the contribution to f\ + fi comes from the region with a s (l/b c ) < 1, or 
equivalently, b c < 0.5/A. At r\ = 0.4, 55% of the contribution is accumulated 
in this perturbative region. As 77 = 1, perturbative contribution has reached 
75%. It implies that the modified PQCD analysis of the decay B — > ttIu in 
the range of 77 > 0.3 is relatively reliable according to the criteria given in [?]]. 
The differential decay rate for the specific case B° — > tt~1 + i> with vanishing 
lepton masses is given by 



d£ 
di] 



l2 G 2 F m 5 r] 3 
768tt 3 



\V ub \ 2 R( V ) = \V t 



ub 



(18) 



Substituting the results of /j into eq. (p!8[), we derive the behavior of R(rj) as 
in fig. 2b, which shows a slow decrease with rj. 

In order to have the full spectrum in 77, we approximate dT/dr] in the 
range of r\ < 0.3 by the soft pion limits of fi [§]]: 



lim R(rj) 



G 2 F m 5 r] f 2 B , 2 
192tt 3 / 2 9bb ** ' 



(19) 
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1.1/ 



B 



is the decay 
15| is the BB*tt coupling con- 



which shows a linear relation with 77. Here fs 
constant of the B* meson, and gBB*n ~ 0.8 
stant. We extrapolate eq. ( |T§D to 77 = 0.3 as shown in fig. 2b, and a good 
match between the soft pion and PQCD predictions is observed. Certainly, 
this extrapolation of the soft pion limit may not be reliable, but the match 
of the two different approaches does justify our calculation to some extent. 

It is then possible to estimate the total decay rate T by combining eq. (PH) 
for 77 < 0.3 with the PQCD predictions for 77 > 0.3. We obtain r k 2.4 x 
10 _11 |\4b| 2 GeV, which corresponds to a branching ratio 0.47 x 10 2 IK,t,| 2 for 
the total width (0.51 ± 0.02) x 10~ 9 MeV of the B° meson |L6 . 
experimental limit on the branching ratio of B° — > 7r~/ + z/ is 3.3 x 10 



'ub\ 

Current 
4 fTl 



We then extract the matrix element \V, 



ub 



2.7 x 10 , close to the value 



0.003 given in the literature [16 
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Figure Captions 



Fig. 1. Graphic representation of eq. (f|). 

Fig. 2. (a) Dependence of /1 + /2 on the cutoff b c for (1) 77 = 0.3, (2) 77 = 0.4, 
and (3) rj = 1.0. (b) Dependence of R{rj) on 77 derived from the modified 
PQCD formalism (solid line) and from the soft pion theorems (dashed line). 
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This figure "figl-l.png" is available in "png" format from: 
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